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Summary 

The research work has advanced the inversion-based guidance theory for: 

• systems with non-hyperbolic internal dynamics; 

• systems with parameter jumps; and 

• systems where a redesign of the output trajectory is desired. 

Non-hyperbolic Internal Dynamics 

A technique to achieve output tracking for nonminimum phase linear systems 
with non-hyperbolic and near non-hyperbolic internal dynamics was devel- 
oped [1]. This approach integrated stable inversion techniques, that achieve 
exact-tracking, with approximation techniques, that modify the internal dy- 
namics to achieve desirable performance. Such modification of the internal 
dynamics was used (a) to remove non-hyperbolicity which is an obstruction 
to applying stable inversion techniques and (b) to reduce large preactuation 
times needed to apply stable inversion for near non-hyperbolic cases. The 
method was applied to an example helicopter hover control problem with 
near non-hyperbolic internal dynamics for illustrating the trade-off between 
exact tracking and reduction of preactuation time. 

Future work will extend these results to guidance of nonlinear nonhyperbolic 
systems. 

Systems with Parameter Jumps 

The exact output tracking problem for systems with parameter jumps was 
considered [2]. Necessary and sufficient conditions were derived for the elim- 
ination of switching-introduced output transient. While previous works had 
studied this problem by developing a regulator that maintains exact track- 
ing through parameter jumps (switches), such techniques are, however, only 
applicable to minimum-phase systems. In contrast, our approach is also ap- 
plicable to nonminimum-phase systems and leads to bounded but possibly 
non-causal solutions. In addition, for the case when the reference trajecto- 
ries are generated by an exosystem, we developed an exact-tracking controller 
which could be written in a feedback form. As in standard regulator theory, 
we also obtained a linear map from the states of the exosystem to the de- 
sired system state, which was defined via a matrix differential equation. The 



constant solution of this differential equation provided asymptotic tracking, 
and coincided with the feedback law used in standard regulator theory. 
Future work will extend these results to nonlinear systems and generalize 
these results to hybrid systems. 

Output Trajectory Redesign 

We studied the optimal redesign of output trajectory for linear invertible sys- 
tems [3, 4]. The specified output trajectory uniquely determines the required 
input and state trajectories, that are found through inversion. These input- 
state trajectories exactly track the desired output, however, they might not 
meet acceptable performance requirements. For example, the required in- 
puts might cause actuator-saturation during an exact tracking maneuver (as 
in the flight control of conventional take-off and landing aircraft). In such 
situations, a compromise is desired between the tracking requirement and 
other goals like reduction of internal vibrations and prevention of actuator 
saturation - i.e., the desired output trajectory needs to be redesigned. 

We posed the trajectory redesign problem as an optimization of a general 
quadratic cost function, and solved it in the context of linear systems. 
Future work will be aimed at extending these results to the nonlinear case 
and at developing algorithms for practical implementation. 
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Abstract - A technique to achieve output tracking for nonminimum phase linear sys- 
tems with non-hyperbolic and near non-hyperbolic internal dynamics is presented. This 
approach integrates stable inversion techniques, that achieve exact-tracking, with approx- 
imation techniques, that modify the internal dynamics to achieve desirable performance. 
Such modification of the internal dynamics is used (a) to remove non-hyperbolicity which is 
an obstruction to applying stable inversion techniques and (b) to reduce large preactuation 
times needed to apply stable inversion for near non-hyperbolic cases. The method is applied 
to an example helicopter hover control problem with near non-hyperbolic internal dynamics 
for illustrating the trade-off between exact tracking and reduction of preactuation time. 

1 Introduction 

Precision output tracking controllers are needed to meet increasingly stringent performance 
requirements in applications like flexible structures, aircraft and air traffic control, robotics, 
and manufacturing systems. While perfect tracking of minimum phase systems is relatively 
easy to achieve, output tracking of nonminimum phase systems tends to be more challenging 
due to fundamental limitations on transient tracking performance [1]. This poor transient 
performance has been mitigated by using preactuation in the stable-inversion based ap- 
proaches for nonminimum phase systems [2, 3, 4, 5]. However, the required preactuation 
time (during which most of the preactuation control effort is required) is large if the zeros of 
a linear system that lie on the open right half of the complex plane are close to the imagi- 
nary axis. In the limiting case, with the zeros on the imaginary axis (non-hyperbolic internal 
dynamics), presently available inversion-based techniques fail because the preactuation time 
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needed becomes infinite. In this paper we present a design technique for output tracking of 
linear non-minimum phase systems, which have non-hyperbolic and near non-hyperbolic non- 
minimum phase internal dynamics. This technique is then applied to an example helicopter 
hover control and simulation results are presented. 

Output tracking has a long history marked by the development of regulator theory for 
linear systems by Francis and Wonham [6] and generalized to the nonlinear case by Byrnes 
and Isidori [7]. These approaches asymptotically track an output from a class of exosystem- 
generated outputs. Further, the Byrnes-Isidori regulator has been extended in [8, 9] and 
computational issues have been studied in [10, 11]. Although nonlinear regulator design is 
computationally difficult, the linear regulator is easily designed by solving a manageable set 
of linear equations. A problem, however, with the regulator approach is that the exosystem 
states are often switched to describe the desired output - this leads to transient tracking- 
errors after the switching instants. Such switching caused transient errors can be avoided by 
using inversion-based approaches to output tracking [4, 12]. Thus, it is advantageous to use 
inversion-based output tracking when precision tracking of a particular output trajectory is 
required. 

Inversion, which is key to our approach, was restricted to causal inverses of minimum 
phase systems in the early works by Silverman and Hirschorn (e.g., [13, 14]) because they 
lead to unbounded inverses in the nonminimum phase case. Di Benedetto and Lucibello [15] 
considered the inversion of time varying nonminimum phase systems with a choice of the 
system’s initial conditions. Instead of choosing initial conditions, preactuation was used by 
noncausal stable inversion techniques developed in [2, 3, 4, 16, 5]. Such noncausal inverses, 
which require preactuation, have been successfully applied to the output tracking of flexible 
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structures [17, 18] and aircraft and air traffic control [19, 20]. However, the fundamental 
limitation of presently available inversion- schemes is that they fail if the internal dynamics 
is non-hyperbolic. Even when the internal dynamics is hyperbolic, if the right-half-plane 
zeros of the system are close to the imaginary axis (the near non-hyperbolic case) then 
the required preactuation time tends to become unacceptably large. In summary, output 
tracking remains a challenge for" nonminimum phase systems with non-hyperbolic or near 
non-hyperbolic internal dynamics. 

There are several approximation based output tracking techniques, where the central 
philosophy is to replace the internal dynamics with a dynamics that provides satisfactory 
behavior, and then to develop the controller based on the altered system [21, 22, 23]. The 
technique most relevant to this paper is developed by Gopalswamy and Hedrick [23] - where 
trajectory modifications are considered to stabilize the internal dynamics. This technique, 
however, requires hyperbolicity of the internal dynamics for computational purposes. The de- 
velopment of computational techniques for stable inversion (e.g., [16]) motivates the present 
integration of the stable inversion scheme with approximation techniques - especially for sys- 
tems with non-hyperbolic internal dynamics where the existing stable-inversion techniques 
fail. However, instead of stabilizing the unstable internal dynamics, we only aim to modify 
the non-hyperbolic behavior with a small perturbation of the internal dynamics. Addition- 
ally, in nonminimum phase systems with near non-hyperbolic internal dynamics the present 
approach allows a tradeoff between the precision tracking requirement and the amount of 
preactuation time, needed to apply the stable inversion based output tracking technique. 

The approximate inversion-based technique is developed in Section 2 and the technique 
is applied to a helicopter hover control example in Section 3, where simulation results are 
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discussed. Our conclusions are in Section 4. 


2 Stable Inversion 

2.1 Inversion-Based Output Tracking Scheme 

Here we describe how the inversion approach is used to develop output tracking controllers. 
Consider a linear system described by 


x(t ) = A x(t) + B u(t ) 


( 1 ) 


y(t) = C x(t) 

which has the same number of inputs as outputs, u(t),y(t) € 3£ p , and x(t) G 3?". Let yd(-) be 
the desired output trajectory to be tracked. Then in the inversion-based approach we, first, 
find a nominal input-state trajectory, [«//(•), x re f(-)] that satisfies the system equations (1) 
and yields the desired output exactly, i.e., 


®re/(^) — A X re j{t^j -|- B Ufj(t) 

Vd{t) = C X ref (^) 

and, second, we stabilize the exact-output yielding state trajectory, x re /( •), by using state 
feedback (see Figure 1). Thus x{t) — > x re f(t) and y(t) — > y<i(t) as t — > oo and output tracking 
is achieved. It is noted that in this output tracking scheme, the reference state trajectory 
x re j ( •) and the feedforward input «//(•) are computed off-line. While stabilization of the 
reference state trajectory can be easily achieved through standard techniques [24] like state 
feedback of the form K[x{t) — x re f{t )], the main challenge is to find the inverse input-state 
trajectory [«//(•), x re f(-)] - especially for systems with nonminimum phase dynamics. This 


| Vi € (-00,00) ( 2 ) 
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paper addresses the off-line computation of the inverse input-state trajectory for a given 
desired trajectory, yd(:)- 


2.2 The Internal Dynamics 

In this subsection, it is shown that finding the inverse input-state trajectory is equivalent to 
finding bounded solutions to the system’s internal dynamics. Let the linear system (1) have 
a well defined vector relative degree, r := [ri, r 2 , . . . , r p ). Then the output’s derivatives are 
given as [25] 


dt rk Vk 


= C k A rk x + C k A rk ' 1 Bu 


(3) 


where C k is the k th row of C, and 1 < k < p. In vector notation let equation (3) be rewritten 


as 


y (r) (i) = A x x(t) + B y u(t) 


(4) 


where 


y (r) := 


d^yi d rz y 2 
dt r * 


d rp y p 

dt r p 


lT 



C x A r ' 


CiA^-'B 


C 2 A r * 


C 2 A r >~ 1 B 

A x := 

C p A rp 

; By 

C v A rp ~ l B 


and By is invertible because of the well-defined relative degree assumption. Equation (4) 
motivates the choice of the control law of the form 


Uff(t) = By 1 [y d (r) (*) - A x x(t)] 


(5) 
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for all t € (— 00,00). Substituting this control law in equation (4), it is seen that exact 
tracking is maintained, i.e., 

y {r] (t) = 

To study the effect of this control law consider a change of coordinates T such that [25] 

m 

»?(*) 


= T x(t) 


'r D -l 


d r ? 


where ((t) consists of the output and its time-derivatives 

d r i - 1 _ d r2_1 

C(0 := [yum, • • • , ’ y*’ * • • ’ dfr 2 — I y2 > • • • » Vp* 3/p» • • • » ■ 

The system equation (1) can then be re-written in the new-coordinates as 

£(*) = AiC + A277 + Bin 

»?(*) = A 3 c + + b 2 u 

where 


(6) 

(7) 


li 

1 

£ 

li 

1 

t— 1 

A2 


A 

A3 


; and B := 

1 

<CQ 

1 


1 



= T-'B 

In the new co-ordinates, the control law for maintaining exact tracking (Equation (5)) can 
be written as 

«(<) = V [»/>(() - AU*) - A<i(t)] (8) 

where 

M 

:= A X T~\ 


Ac 

A„ 


Note that the desired £(•) is known when the desired output trajectory yd(-) and the output’s 
time derivatives are specified. This desired ((■) is defined as Q(-). Since the control law was 
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chosen such that exact tracking is maintained, y^(t) = y^(f) we also have £(i) = and 

Equations (6) and (7) become 


m 

= C At) 

(9) 

v(t) 

= A 3 c + A t ri + B 2 B [»/>(() - A c Q(t) - AflM] ■ 

(10) 


This is the inverse system, and in particular, Equation (10) is the internal dynamics. Solving 
the internal dynamics is key to finding the inverse input-state trajectories. If a bounded 
solution, rid(’), to the internal dynamics (10) can be found then the feedforward input can 
be found through equation (8) as 


“//(*) = B y 1 [2/d (r) (*) ~ ArCKO “ Ar,r} d {t)] (11) 

and the reference trajectory can be found as 


x ref (t) = T 1 


Cd(t) 


[ J 

Thus a bounded solution to the internal dynamics (10) is required to find the inverse and to 
apply the output tracking scheme shown in Figure 1. 


2.3 Modified Internal Dynamics 

Standard inversion schemes [13, 14] that integrate (forward in time) the internal dynam- 
ics (10) lead to unbounded solutions since the internal dynamics is unstable for nonmini- 
mum phase systems. Noncausal inversion (e.g., [4]) leads to bounded but noncausal solution 
to the internal dynamics. Such stable-inversion techniques are, however, not applicable to 
systems with non-hyperbolic internal dynamics. In this subsection a compromise between 
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stable inversion and approximation based inversion schemes is proposed. The key is to 
modify the internal dynamics by giving up exact output tracking - enough to remove the 
non-hyperbolicity, and then to apply stable-inversion. The difference between the proposed 
technique and other approximation techniques is that the internal dynamics is perturbed 
only to remove the non-hyperbolicity, and not to stabilize the entire internal dynamics. 

To change the non-hyperbolicity of the internal dynamics an extra term, v(t), is added 
to the control law (8) as follows 


u ff (t) = B y 1 [y/ r \t) - A c ( d (t) - A v r/(t) + u(t)] . 


( 12 ) 


With this modified control law the inverse system (Equations (9) and (10)) become 


d_ 

dt 


e c (t) 

= s 

1 

■'o ’ 1 

CD 

1 

1 

CS- 

1 


. *>(*) . 


+ G y Yd(t) + G v v(t ) 


(13) 


where e^(t) := £(t) — Crf(^) is the error in the output and the output’s derivatives, 

0 0 
i 3 (a, + B 2 B^A V ) J 


5 


J Gy = 


BiB - 1 


BoB~ l 


Gy 


, and Y d (t) 


Cd(t) 


0 0 
[ (a - fyB^Ac) B 2 B - 1 J 

Assuming that the original system (A, B) is controllable, we also have (S,G V ) controllable 
and hence a there exists a feedback of the form 

e c {t) 

n(t) 


v{t) = F 


(14) 


such that the modified inverse system (13) is hyperbolic — i.e., all poles on the imaginary axis 
are moved. Note that this change to an hyperbolic system can be achieved through arbitrarily 
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small F since non-hyperbolicity is not structurally stable property. The hyperbolic system 


d_ 

dt 


e c (t) 

v(t) 


(S + G V F) 


e c {t) 

*7(<) 


+ GyYd(t) 



+ GyY d (t) 


(15) 


is the modified inverse system. This modification of the internal dynamics can also be 
used to move unstable poles of the inverse system that may be close to the imaginary axis 
for reducing the required preactuation time. Next, stable inversion of the modified inverse 
system is carried out [4]. 


2.4 Computation of the Inverse 

We begin by decoupling the modified internal dynamics (15) into stable (z s ) and unstable (z u ) 
subsystems. Since the modified internal dynamics is hyperbolic, there exits a decoupling 
transformation U such that the modified inverse system (15) can be written as 


z s {t) = S s z s (t) + G s Yd(t) 
z u {t) = S u z u (t) + G u Y d (t) 


(16) 


where 


z(t) := 

z s(t) 


t 

i 

= u 


z u(t) 


v(t) 


(17) 


To find bounded solutions to the unstable inverse systems, the boundary conditions that 
z s (— oo) = 0 and z u ( oo) = 0 are applied to Equation (16). This leads to unique bounded 
solutions to the modified internal dynamics by flowing the stable subsystem forward in time 
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and flowing the unstable system backward in time - as 


z s ,d{t) = Jloo e s ^ t ~ T '>G u Y d {T)dT Vie (- 00 , 00 ), 


( 18 ) 


z u4 {t) = - f t °° e s ^- T )G u Y d (r)dr Vi € (- 00 , 00 ). 

This completes the technique. To summarize, the bounded solution (18) is used to find the 
reference state trajectory as x re} (t) - T~ l U~ l z d {t), and to find the feedforward input, 
Ufj(-) from equation (12). This inverse, [ujf(-), z re /(-)], is then used in the control scheme 
shown in Figure 1 to obtain output tracking. 


2.5 Preactuation Time and Unstable Poles of Inverse System 

The connection between the amount of preactuation time required to apply the inversion- 
based feedforward input and the unstable poles of the modified inverse is established in the 
following Lemma. Lemma Let 

1. the support of Y d (-) be contained in [to, 00 ) for some to, 

2. all the unstable poles of internal dynamics represented by the eigenvalues of S u lie to 
the right, in the complex plane, of the line Re(s) = a for some positive a, and 


3- II G.Y4-) IL < /?• 


Then there exists M such that | \ufj(t) < Me a( - t ~ to ' > for all time before the start of the 
maneuver, t < to- 

Proof: From Condition 2 of the Lemma, there exists a positive constant Ms u such that 


e Su(t-r) < M Su e a ^ Vi < t . 


(19) 
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Then for all t < i 


|| u ff (t) IL = \\B; 1 [y d {r \t) - A c Cd(t) - A n i)(t) + 1>(*)] L from Equation (12) 

= l^" 1 [— A v r\(t) + u(t)] | from Condition 1 of the Lemma 

1 \\ 1 

0 0 e c (t) 

+ F from Equation (14) 

0 -B^Ar, J T7(t) 

J f L. J OO 

e c (t) 

Au 

v{t) 

J CO 

e c (t) 

V (t) 

J OO 

Z s (t) 

< H^IUICf- 1 U from Equation (17) 

Zu(t) 

J OO 

= II Au II ^ \\U~ l HJM*) IL Since z a (t) = 0 for all t < t 0 
= II Au ILII^- 1 |L||r e s “f < - T )G u ld(r)dr ^ from Equation (18) 

= II Mi IUIC7 -1 lUll/- e s “ (<_T) G u Fd(r)dr ^ from Condition 1 of Lemma 

< PM Sv \\Au IUII7” 1 IL r e“ (t - T) )dr 

J to 

from Equation (19) and Condition 3 of Lemma. Integrating the above expression we get 

IK/(<) IL < lUltr - 1 

:= Me a(t ~ to) , 





which concludes the proof. ® 

The above Lemma states that the preactuation input tends to zero exponentially, as we 
go back in time from the start of the maneuver at to. The rate at which the preactation 
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becomes zero can be increased by moving the unstable poles of the modified inverse away 
from the imaginary axis - at the expense of exact output tracking. The trade-off between 
exact-tracking of the desired output and reduction of the pre-actuation time is illustrated in 
the following example. 


3 Example: Helicopter Hover Control 


Here, we apply the output-tracking technique to the hover control of a Bell 205 Helicopter, 
which has near non-hyperbolic unstable internal dynamics. We consider one of the cases 
studied previously in [26] , wherein the aircraft dynamics was trimmed at a nominal 5° pitch 
attitude, with a mid-range weight, a mid-position center of gravity and operating in-ground 
effect at near sea level. The linearized model is given as [26, 27] 


x 


where 


0 

0.03 

0.18 

-0.10 

-0.39 

0.09 

0.01 

-0.01 

-0.19 

0.02 

0 

-0.41 

0.03 

-0.02 

-0.88 

-0.01 

-0.02 

-0.06 

0 

0 

1 

0 

0 

0 


= Ax + Bu 


-0.01 

-0.42 

0.08 

-9.81 

0 

-0.10 

-0.72 

0.68 

0 

0 

0 

0.23 

0.04 

0 

0 

-0.05 

-0.27 

0.27 

0 

9.81 

-0.04 

-0.57 

0.14 

0 

0 

0.07 

-0.32 

-0.71 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 


( 20 ) 


( 21 ) 
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0.08 0.13 


0 


0 


- 1.17 

0.04 

0 

0.01 

0 

- 0.07 

0 

0.01 

- 0.04 

0 

0.11 

0.19 

- 0.04 

0 

0.22 

0.17 

0.17 

0 

0.03 

- 0.47 

0 

0 

0 

0 

0 

0 

0 

0 


U 

W 

Q 

V 

p 


forward velocity 
vertical velocity 
pitch rate 
lateral velocity 
roll rate ) 


R 

6 


yaw rate 
roll rate 
pitch attitude 


* 


L 

J 

roll attitude j 



collective 

&B 


longitudinal cyclic 

8a 


lateral cyclic 

Sp 


tail rotor collective 


( 22 ) 


( 23 ) 


( 24 ) 



It is noted that the helicopter’s actual dynamic behavior differs because of modeling er- 
rors like nonlinearities and unmodeled dynamics. In output tracking control schemes that 
depend on the model, such modeling errors need to be corrected through feedback in the 
control scheme (see Figure 1). In particular, modeling errors can be compensated by robust 
stabilization of the reference state-trajectory, see for example [26]. The goal of the present 
paper is to develop inversion-based feedforward and reference state trajectories for use in 
the control scheme shown in Figure 1. In the following, we apply the inversion technique to 
control the helicopters forward, lateral and vertical velocities and its yaw rate. The forward 
velocity and the yaw rate are to be kept at zero and the desired profiles of lateral and vertical 
velocities and accelerations are as shown in Figures 2 and 3. 

3.1 Internal Dynamics 

To find the internal dynamics we begin with a change in the co-ordinates. Let, ( be defined 
as the outputs 

m 

W(t ) 

V(t) 

m 

and let r) be the remaining states, 

Q(t) 

m 

(26) 

p(t) 

4>{t) 





14 



In the (£, 77) co-ordinate system, given by 


1 0 0 0 0 0 0 0 

0 1 0 0 0 0 0 0 

0 0 0 1 0 0 0 0 

C(t) . 0 0 0 0 0 1 0 0 

x(t ) 

rj(t) J 0 0 1 0 0 0 0 0 

0 0 0 0 0 0 1 0 

0 0 0 0 1 0 0 0 

0 0 0 0 0 0 0 1 

:= Tx(t) , (27) 

the system equations can be re-written as 

C(t) 

= TAT - 1 + TBu(t) 

. . 

Ai A 2 ((t) B x 

+ u(t) . (28) 

A 3 A 4 rj(t ) B 2 

Given a desired output trajectory and the desired output’s time derivatives, [CdO)? CdOL the 
exact output tracking control law (see Section 2.4) is found from equation (28) as 

Uff(t) = Bi 1 [OCO - AiCd(t ) - A 2 T](t)\ (29) 

With this control law the inverse system becomes (from equation (28)) 

C(t) = (d(t) ( 30 ) 

r)(t) = A 4 r/(t) +A3C(t) + B 2 u(t) 
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= [a, - b 2 b~ 1 a 2 \ rj(t) + B 2 B~ l [6(0 - AMt)] 
:= A v T](t) + B 2 B 1 1 [Crf(^) ~ AC<A)1 • 


( 31 ) 


The problem is solved by finding a bounded solution to the internal dynamics (31). However, 
the bounded solution found through stable-inversion is noncausal and could require a large 
preactuation time if the poles of the internal dynamics are unstable and lie close to the 
imaginary axis in the complex plane. For this particular example, there are two such complex 
conjugate poles near the imaginary axis, 0.0425 db 4.3055*. We modify the exact-tracking 
control law (29) to shift these poles away from the imaginary axis to 2 ± 4.3055*. This is 
described next. 

3.2 Modified Inverse System 

Following the approach described in Sections 2.3, we modify the internal dynamics by adding 
a term v(t ) to the control law (29) to obtain 

u ff (t) = B- 1 [6(t) - AMt) - A 2Vd {t)] + v(t) . (32) 

Substituting this control law into Equations 30 and 31, the modified inverse system is ob- 
tained as 

e^(t) 0 0 e^(t) Bi 

= -1- v(t) + 

V(t) A Ar, rj(t) B 2 

0 0 C lit) 

(A - nr 1 A) B 2 Bx 1 J [ Ut) 
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A e c (t) 

:= 5 + G v v(t) + G y Y d (t). (33) 

. . 

where e^(t) := ((t) — (d{t)- The poles of the inverse system can be moved to any desired 
location by using the control v(t) because (S', G v ) is controllable. However, such modifica- 
tions, aimed at reducing preactuation time, will also lead to a loss of precision in output 
tracking. This tradeoff between the reduction of preactuation time and the loss of precision 
in tracking is illustrated through simulation. 

3.3 Simulation Results and Discussion 

Two sets of simulations were performed. First, stable inversion was applied to the orig- 
inal system which lead to exact output-tracking inverse input-state trajectories. Second, 
simulations were performed when the unstable poles of the inverse system are moved from 
0.0425 ± 4.3055? to 2 ± 4.3055? for reducing the amount of preactuation time required. 
Further, the inverse system also has four poles at the origin - corresponding to four pure 
integrators for (e<;(-)) dynamics - which were moved to — 1,-2, —3 and —4 for stability of 
the numerical integration scheme. 

Figures 2 and 3 show the desired output trajectories for the lateral and vertical motions 
(corresponding to unit displacements in the two directions), while the forward velocity and 
yaw rate were to be maintained at zero - the maneuver starts at time t = 0. Figures 4 and 5 
show the output trajectories achieved by the inverse state-trajectory, x re /(-), which is to be 
used as a reference trajectory in the feedback scheme shown in Figure 1. The corresponding 
feedforward inputs are shown in figures 6 and 7. Note here that the feedforward inputs 
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are non-zero before the start of the maneuver, i.e., time t < 0, and hence preactuation is 
required. 

Figures 4 and 5 show that exact-output tracking reference state trajectories can be found, 
even when the internal dynamics is unstable, through the stable inversion approach. The 
stable-inversion technique yields bounded solutions to the unstable internal dynamics, i.e., 
the pitch and roll motions are bounded as shown in Figure 8. However, the feedforward input 
found through exact-inversion requires substantial preactuation time as shown in Figures 6 
and 7, i.e., the preactuation remains non-zero for a significant time before the start of the 
maneuver at t = 0. Figure 9 shows that about 30 seconds of preactuation is needed to 
apply the inverse of the original system for output tracking - in contrast, modification of the 
internal dynamics reduces the preactuation needed from 30 seconds to 1 second (see Figure 
9). As seen in Figures 4 and 5 the output trajectories are still tracked well by the modified 
inverse. Further, this substantial reduction in preactuation time is achieved with similar 
control efforts and with similar roll and pitch motions (see Figures 6-8). Thus the approach 
presented here, allows a trade-off between precision tracking and the amount of preactuation 
that is acceptable. Future work will generalize the results to nonlinear nonminimum phase 
systems with non-hyperbolic internal dynamics. 

4 Conclusions 

A technique to achieve output tracking for nonminimum phase linear systems with non- 
hyperbolic and near non-hyperbolic internal dynamics was presented. This approach is an 
integration of the stable inversion technique that aims at exact tracking with the approxi- 
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mation approach that modifies the internal dynamics to achieve desirable performance. The 
method was applied to an example helicopter hover control problem to illustrate the effects 
of modifying the internal dynamics. It was shown that substantial reduction in preactuation 
time is possible by giving up some of the precision in tracking - thus making the stable 
inversion approach viable for practical application. 

References 

[1] L. Qui, and E. J. Davison. Performance limitaions of non-minimum phase systems in 
the servomechanism problem. Automatica, 29, March:337-349, 1993. 

[2] E. Bayo. A finite-element approach to control the end-point motion of single-link flexible 
robot. J. of Robotic Systems, 4,(l):63-75, 1987. 

[3] B. Paden and D. Chen. A state-space condition for the invertibility of nonlinear 
nonminimum-phase systems. Advances in Robust and Nonlinear Control Systems, DSC- 
ASME, 43:37-41, 1992. 

[4] S. Devasia, D. Chen, and B. Paden. Nonlinear inversion-based output tracking. IEEE 
Transactions on Automatic Control , 41(7):930-943, 1996. 

[5] L.R. Hunt, G. Meyer, and R. Su. Noncausal inverses for linear systems. IEEE Trans- 
actions on Automatic Control , 41 (4) :608— 61 1 , 1996. 

[6] B. A. Francis and W. M. Wonham. The internal model principle of control theory. Au- 
tomatica, 12:457-465, 1976. 


19 



[7] A. Isidori and C. I. Byrnes. Output regulation of nonlinear systems. IEEE Transactions 
on Automatic Control, 35(2):131— 140, 1990. 

[8] B. Castillo. Output regulation of nonlinear systems with more inputs than outputs. 
International Journal of Control, 57,(6) :1343— 1356, 1993. 

[9] J. Huang. Output regulation of nonlinear systems with nonhyperbolic zero dynamics. 
IEEE Transactions on Automatic Control, 40(8) :1497— 1500, 1995. 

[10] J. Huang and W. J. Rugh. An approximation method for the nonlinear servomechanism 
problem. IEEE Transactions on Automatic Control, Sept:1395-1398, 1992. 

[11] A. J. Krener. The Construction of Optimal Linear and Nonlinear Regulators. In Sys- 
tems, Models and Feedback: Theory and Application. A. Isidori and T.J. Tam (editors). 
Birkhauser, Boston, 1992. 

[12] S. Devasia, B. Paden, and C. Rossi. Exact-output tracking theory for systems with 
parameter jumps. International Journal of Control, Accepted, 1996. 

[13] L. M. Silverman. Inversion of multivariable linear systems. IEEE Transactions on 
Automatic Control, 14(3):270-276, 1969. 

[14] R. M. Hirschorn. Invertibility of multivariable nonlinear control system. IEEE Trans- 
actions on Automatic Control, 24(6):855-865, 1979. 

[15] M. D. Di Benedetto and P. Lucibello. Inversion of nonlinear time- varying systems. IEEE 
Trans. Automatic Control, 38:1259-1264, 1993. 


20 



[16] S. Devasia and B. Paden. Exact output tracking for nonlinear time-varying systems. 
Proceedings of the 33rd CDC, Lake Buena Vista, Florida , pages 2346-2355, 1994. 

[17] B. Paden, D. Chen, R. Ledesma, and E. Bayo. Exponentially stable tracking control for 
multi-joint flexible manipulators. ASME Journal of Dynamic Systems, Measurement 
and Control, 115:53-59, 1993. 

[18] S. Devasia and E. Bayo. Redundant actuators to achieve minimal vibration trajec- 
tory tracking of flexible multibodies: Theory and application. Journal of Nonlinear 
Dynamics , 6:419-431, 1994. 

[19] P. Martin, S. Devasia, and B. Paden. A different look at output tracking: Control of a 
vtol aircraft. Automatica , 32(1):101-107, 1996. 

[20] G. Meyer, L. R. Hunt, and R. Su. Nonlinear system guidance in the presence of trans- 
mission zero dynamics. NASA Technical Memorandum 4661, January, 1995. 

[21] R. Gurumoorthy and S. R. Sanders. Controling non-minimum phase nonlinear sys- 
tems - the inverted pendulum on a cart example. Proceedings of the American Control 
Conference, pages 680-685, 1993. 

[22] Z. Xu and J. Hauser, ’’Higher order approximate feedback linearization about a manifold 
for multi-input systems”a, IEEE Transactions on Automatic Control. Vol. 40 No.5, pp. 
833-840, 1995. 

[23] S. Gopalswamy and J.K. Hedrick. Tracking nonlinear non-minimum phase systems 
using sliding control. International Journal of Control, 57.5:1141-1158, 1993. 


21 



[24] Hassan K. Khalil. Nonlinear Systems. Macmillan, 1991. 


[25] A. Isidori. Nonlinear Control Systems: An Introduction. Springier- Verlag, 1989. 

[26] J. K. Pieper, S. Baillie, and K. R. Goheen. Linear quadratic optimal model-following 
control of a helicopter in hover. Optimal Control Applications and Methods , 17:123-140, 
1996. 

[27] R. K. Heffley, W. F. Jewell, J. M. Lehman, and R. A. Van Winkel. A compilation of 
helicopter handling qualities data. NASA Contract Report 3144 under Grant NAS2- 

9344, K 1979. 


22 




Figure 1. The Control Scheme 
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Figure 2. Desired Lateral and Vertical Velocity Profile. 
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Figure 3. Desired Lateral and Vertical Acceleration Profile. 
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Figure 4. Lateral and Vertical Velocity achieved by the Inverse Reference State Trajectory. 
The dotted line is for the exact-tracking case without modification of the internal dynamics 

and the dashed line is with modification. 
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Figure 5. Forward Velocity and Yaw Rate achieved by the Inverse Reference State 
Trajectory - the desired value is zero. The dotted line is for the exact-tracking case 
without modification of the internal dynamics and the dashed line is with modification. 
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tail rotor collective lateral cyclic 




Figure 7. Feedforward Inputs. Dotted line is without modification of internal dynamics 

and the dashed line is with modification. 
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Figure 9. Comparison of Required Preactuation in the Feedforward. 
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Abstract - In this paper we consider the exact output tracking problem 
for systems with parameter jumps. Necessary and sufficient conditions are 
derived for the elimination of switching-introduced output transient. Previ- 
ous works have studied this problem by developing a regulator that main- 
tains exact tracking through parameter jumps (switches). Such techniques 
are, however, only applicable to minimum-phase systems. In contrast, our 
approach is applicable to nonminimum-phase systems and obtains bounded 
but possibly non-causal solutions. If the reference trajectories are generated 
by an exosystem, then we develop an exact-tracking controller in a feedback 
form. As in standard regulator theory, we obtain a linear map from the 
states of the exosystem to the desired system state which is defined via a 
matrix differential equation. The constant solution of this differential equa- 
tion provides asymptotic tracking, and coincides with the feedback law used 
in standard regulator theory. The obtained results are applied to a simple 
flexible manipulator with jumps in the pay-load mass. 

1. Introduction 

We study the exact-output tracking of systems, which are described by 

f *(*) = A[fc(t)]a:(t) + £[fc(<)]u(t) 

\ y{t ) = C[k(t))x(t) ' > 

where x € 3? n , with same number of inputs as outputs u(t),y{t) € 9? p . The 
system matrices A(k), B(k),C(k ) are constant over time intervals /*, where 
k belongs to a finite index set tC ^ [0, and the parameter change 
(switch) at times t = (See Figure 1). Here, the switching time 

are known, in contrast to systems where the switches may be signal driven. 

For constant linear systems, the asymptotic output-tracking problems for 
constant linear systems has received much attention in the past. In particu- 
lar, the regulator theory ( Francis (1977); Basile and Marro (1992); Wonharn 
(1985)) provides a general framework in which the asymptotic output track- 
ing can be solved when the reference trajectory is generated through a linear 
exosystem. In the presence of switches in the system, one technique for 
achieving output regulation is to switch the regulator. Note that regulation 
can be recovered between two consecutive switches (due to asymptotic prop- 
erties) - especially if the switching occurs far apart in time. However, this 
technique also tends to induce transients in the output during the switches. 
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Figure 1: The Switching Times 


In order to eliminate these switching-caused transients, a regulation scheme 
that maintains exact trajectory tracking across system switches must be used. 
This fairly new problem has been studied in ( Marro and Piazzi (1993)), for 
minimum-phase systems. In this work, a feedforward action is used in con- 
junction with the feedback defined by the regulator in order to cancel the 
output transients across the switches. 

In the present paper we propose an alternative approach for exact output- 
tracking of switched-systems, which is also applicable to nonminimum-phase 
systems. In the nonminimum-phase case, a bounded non-causal solution 
is obtained ( Devasia et al. (1996)) that requires the preknowledge of the 
reference trajectory and of all the switching times. Such exact tracking 
schemes based on noncausal schemes is useful in problems like aircraft guid- 
ance ( Meyer et al. (1995) and Hunt et al. (1996)). 

We present necessary and sufficient conditions for the solvability of the in- 
version problem for linear systems with switches; the inverse is used to track 
the desired output. We consider two kinds of desired output trajectories: 

(a) a single pre-specified trajectory or (b) belonging to a class of outputs gen- 
erated by a given linear exosystem, that could undergo parameter changes 
as well. In this latter case, we obtain the solution in a time-varying feed- 
back form, where the feedback matrix satisfies a matrix ordinary differential 
equation. The equilibrium solution of this differential equation solves the 
asymptotic output tracking problem, and coincides with the feedback ma- 
trix resulting from the standard regulator. This establishes an interesting 
connection between our approach and the traditional regulator theory. 

The paper is organized as follows: in Section 2 the exact tracking of a 
single prescribed output trajectory is considered, and necessary and sufficient 
conditions are presented. A geometric version of the obtained conditions is 
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also provided. In Section 3, the case of reference trajectory obtained through 
a linear exosystem will be treated. The conditions of the previous section 
when rearranged establish a close relationship with the traditional theory of 
output regulation. Section 4 focuses on the additional problem of stability 
of the closed loop system. Finally, section 5 presents the application of the 
developed theory to a simple nonminimum-phase switched system, given by 
a flexible beam subjected to step variation of the pay-load mass. Conclusions 
end the paper. 

2. Tracking a prescribed output trajectory 

Below we formulate the exact tracking problem for a prescribed output tra- 
jectory, and establish necessary and sufficient conditions for its solvability. 
Geometric interpretations of these conditions are also provided. 

The Inversion Problem 

Given a desired output trajectory yd, find a pair of state and input tra- 
jectories Xd, Ud such that 

1. Xd,Ud satisfy the system equation (1): 

= ^4[fc(<)]x,i(<) + B[k(t)]ud(t); Vt G (— oo, oo), (2) 

2. exact output tracking is achieved (even across switches): 

Vd(t) = C'[fc(<)]x«j(<); Vi € ( oo, oo), 

3. and the inputs and state are bounded: 

II *<*(•) lloo < 00 > 

II **<•) lloo < 00 • 

Using the Inverse for Exact- Output Tracking 

The existence of an inverse (ud,Xd) implies that there are input-state 
trajectories that yield the desired output - exact output tracking is easily 
achieved by stabilizing the desired state trajectory. For example, use Ud as 
feedforward and use the error x — Xd for feedback (see Figure 2). Stabilization 
is not the central issue in this paper, and any scheme for feedback design can 
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Figure 2: The Control Scheme 

be used. For example, given (A(k), B(k)) controllable for all k, the system 
may be stabilized through pole placement with all the closed loop poles in 
the same locations for all k. 

Note that typically, the initial conditions of the system are different from 
the initial conditions of the desired state trajectory leading to initial transient 
errors typical of all tracking controllers. However, once the desired level of 
tracking is achieved (due to exponential reduction in error) our technique will 
maintain tracking across parameter switches. In contrast, switching standard 
regulators when the system parameters change will cause transient errors at 
the switching instants - tracking will not be maintained across switches. 


Exact- Tracking Maintaining Input 
We will assume the following: 


Assumption 1 System A(k), B(k),C(k) has well defined vector relative de- 
gree ( Isidori (1989)) for each k € 1C. 

Then we can find a co-ordinate transformation Qk such that ( Isidori (1989)) 


Y k (t) 

Zk(t) 


= Q k x{t) = 


c k 

Z k 


x{t) 


where 

%(*) = [yi, yi, • • • , [j yi» y2, m, 


d(n, 2 -i) dKp- 1 ) 

s/2, • • • , y P , y P , • • • > rfi (r fc>p -i) Vp\ ’ 


and rjt = [rjt,i , rjt >2 , . . . , r^p] is the systems vector relative degree for fjt < 
t < tk+i (where we assume to = — oo). Here, maps the system states into 
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the outputs and its time derivatives and Z k maps the states into the internal 
dynamics, z k . 

Note that a necessary condition for exact output tracking in the interval 
Ik is that the system state at time tk is such that 

n(**) = n,«*(**) 

In addition, to maintain exact tracking we need to ensure that 


(r k ) r dM dM dM y 

yd ~ d<( r *,2) yd ’ 2 ’ d& k ^ yd,p ’ 

then we can find the following unique control law (provided assumption 1 is 
satisfied) ( Isidori (1989)) 

u d {t) = F k x{t) + G k y ( J k) {t) (3) 

such that the time derivatives of the output is the same as that of the de- 
sired output trajectory yd (this is also a necessary condition for exact output 
tracking). This exact-tracking control is completely determined by the state 
x(t k ), and by the desired output along with its derivatives up to the order 

r k - 

Substituting control law (3) in (2) we obtain for t € I k : 

x{t) = A y (k)x(t) + B y (k)y% k \t) (4) 

where i4 y (A;) = A(k ) + B(k)F k , B y {k ) = B(k)G k . In the transformed co- 
ordinates the system equations for t € I k are of the form 

nw = km 

h(t) = A,(k)z t (i) + A,,y(k)Y k (t) + B z (k)y f ;-\t) [ 1 

Our objective is to define under which conditions it is possible to define a 
feasible state trajectory Xd(t) such that exact trajectory tracking is preserved 
through all the time intervals. There are two main hurdles. Firstly, the 
existence of at least a state trajectory which maintains exact output tracking 
has to be determined. This depends on compatibility of the desired output 
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with the given system. Secondly, the state trajectories have to be bounded. 
In systems with unstable internal-dynamics (nonminimum phase systems) 
generic solutions, tend to be unbounded. In this case, we need to establish 
additional conditions for the existence of bounded solutions to the internal 
dynamics. 

In this paper, we restrict ourselves to the case where (a) Yk,d, the output 
along with its time derivatives, has a compact support [T,-, Tj] € (—00, 00), 
(b) the switching occurs within this compact set, and (c) the internal-dynamics 
is hyperbolic. More formally 


Assumption 2 The desired output trajectory yd(') and its time derivatives 
are bounded and have compact support S{ya) = [T,-, T/] . The switching in the 
system parameters are at fixed times tk € (T,-, Tj ) for every k 6 [1, . . . , N}. 

Assumption 3 System (1) has hyperbolic internal dynamics, i.e. the eigen- 
values of A z (k) have non-zero real parts (no centers) for k = 0 and k = N. 
This is equivalent to requiring that the original system (1) have no zeros 
which lie on the imaginary axis ( Isidori (1989)) for k = 0 and k = N. 


The last assumption implies the existence of transformations Qo,Qn such 
that system state can be partitioned into 


’ Y k (t) ■ 


Cl 1 

Zskifi) 

= Q k x(t ) = 

Z a k 

. Z uk(t) _ 


. ^uk . 


( 6 ) 


where z s k and z u k are the coordinates for the stable and the unstable sub- 
spaces of the system’s internal-dynamics. 

Notations 

Towards establishing conditions for the existence of solutions to the ex- 
act tracking problem for a prescribed output, we first study the dynamic 
evolution of the system for a given initial condition. 

Given an initial condition in an interval Ik the system’s evolution for 
Ik < t < tk+i is described by 

Xd(t) = :d(tk) + [ B y {k)y^ h \r)dT . 

Jtk 
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In a more compact form 


Xd(t) = $k(t,tk)xd(tk) + h k (t,t k ) 

where 

$k(M*) = e Av ^*~ tk) 

and 

Jt k 

The above equations describe the flow in an interval where the system 
doesn’t undergo switches. In order to obtain a representation of the system 
state in terms of an initial state that does not belong to the same interval, 
we define flow compositions as follows: 

W k ,i(t,ti) = $ k (t,t k ) 0 $ k -i(t k ,t k -i) 0...0 

k - 1 

Hk,i(t,U) - h k (t,t k ) + ,j(tk-utj)hj(tj +u tj) 

i=* 

where 

and 

U) = hi(t,U) 

The system evolution for an initial condition x(T{) can be rewritten as: 

x(t) = **,o (t, Ti)x(Ti) + H k ,o(t, Ti ) (7) 

Necessary and Sufficient Conditions 
We first formally state the result. 

Lemma 1 Under assumptions 1-3, the exact output tracking problem is solv- 
able with bounded solution if and only if the system of equations: 

[ 0 

Ydk(tk) = ClQkVk-i,o{tk,Ti)Qn l 0 + 

z u,o{Ti) 

C k Q k H k -i )0 (t k ,Ti ) Vk = l,...,N (8) 
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0 = Z uk Q N ^Ar,o(I/,r,)Qo 1 

admits solution in z Ul o(Ti). 

Proof. System trajectories outside [T,-,Ty], the compact support of Vj, are 
bounded if and only if z u n, the unstable component of the internal-dynamics, 
is zero at the end of the motion Tj and similarly the stable component z„o is 
zero before time T,. Formally: 


0 

0 


+ H N ,o(T f , Ti) 


(9) 


•(Ti) = Qo 1 


0 

0 




*(T,) = Qn 1 


0 


z *,n(T/) 

0 


Substituting the previous expressions in (7) computed at t = Tj 

x(Tj) = * Nfi (T f ,Ti)x(Ti) + HsAT/iTi) 

gives (9). In addition, exact tracking in every interval I k is possible if it is 
possible to find state trajectories that are continuous and such that C k x(t k ) = 
Yd k (t k ). By using again (7) that gives the state at t = t k as a function of the 
initial state and the constraint on x(ZJ), (8) easily follows. ■ 

Equation (8) will be referred in the following as compatibility conditions , 
while equation (9) will be referred to as stability condition. The compatibility 
conditions ensure that Yd does not jump across switches (else unbounded 
inputs will be required)! The stability condition ensures that the autonomous 
system dynamics for t — > ±oo is bounded. 

The algebraic conditions expressed by lemma 1 can also be interpreted in 
a geometric coordinate-free framework. To this end, let 


C k = {x :Y k = Yd, k (t k )} k € [0, . . . , N] 


represents the set of the admissible system states at time t k in order to 
achieve exact tracking in the time interval t k < t < t k +i. C k is in general 
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a linear variety in the state space, that reduces to a linear subspace (the 
system internal dynamics) for k = 0. 

A necessary condition for achieving exact tracking when t < T,- is x e j(T t ) € 
£ 0 . Further, to maintain a bounded solution for all t < T, it is necessary 
that the initial state belong to the unstable subspace of the system internal 
dynamics Xd(Ti ) € £«,<>. 

Note that every Xd(T{) determines a unique Xd{ti), given by 

Xd{h) = $o(^i ,Ti)xd(Ti) + ho(ti,Ti) 

Hence, we can define the image of the subspace £ U) o as: 

$o(ti,to) o Cu,o = {x : x = $o{ti,to)y + ho(ti,to)] V € £u,o} 

which represents the linear variety composed by the points reachable at time 
ti with the constraint of y(t) = yd(t ) for all t € [to^i)- 

To maintain exact tracking in the next interval t € 1 , ^ 2 ) ? it is necessary 

that Xd(t 1 ) € C\. The compatibility condition at time t — t\ states that 

Xd(ti) € Ci PI * 0 (*i, 2<) 0 £u,a 

which is possible if and only if the linear variety: 

Si = Ci H $o(^i ,Ti) o C Ut o 

is not empty, i.e. if and only if C\ intersects the image of £ U)0 under the 
system flow. The same procedure can be repeated for the switching time 
t = t 2 . Starting from Si, we can flow forward in time. In order to achieve 
exact tracking in the interval [< 2 ,< 3 ), the image of Si must intersect C 2 , i.e. 
the set 

s 2 = £ 2 n ^i(t 2 >^i) ° <^i 

must be not empty and more generally, the exact tracking in t £ [T,-, tk+i ) is 
possible if and only if the image of Sk-i under the system flow tk-i) 

intersects £*., i.e. the set: 

Sk = Ck n o Sjt_i (10) 

is not empty for every k = 1, . . . , N. However, in order to obtain a bounded 
solution for t > T/, the final state at time t = Tj must belong to the stable 
subspace of the system internal dynamics £ s ,jv. Let 

Sr, A £ s ,w n $tv(T/, tjv) o Sn- 

Hence, we have proved an analogue of lemma 1 in geometric terms 
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Lemma 2 The exact output tracking problem is solvable if and only if St, 
is non-empty, i.e., 

St, * 0. (11) 

Since the flow of a linear system is a homeomorphism, the dimension of 
a linear variety and that of its image are equal, and hence: 

dim(5,) > dim(«Sj); j > i 

This means that at each iteration (10) the set of possible solutions could 
reduce at each K and that no solution is possible if it becomes the empty set 
for some k - i.e., it is empty for every j > k. 

Switched systems with invariant internal-dynamics subspace 

We present below a particular case in which the given conditions consid- 
erably simplify. This exemplifies the obtained results and will be illustrated 
with an example in section 5. Let 

Assumption 4 System (1) has constant relative degree r = r* for every k 
and matrix C% — C* is constant for every k. 

It follows from the previous assumption that the coordinates outside the 
internal dynamics Yjt are the same for every k, and thus the internal dynamics 
subspace is the same for every k. Note that the stable and unstable subspaces 
may still be different, and may switch around, but are constrained to belong 
to the same subspace! 

Since Yk = C*x, the continuity of x implies that the compatibility condi- 
tions are always satisfied. Formally. 

Lemma 3 If assumption 4 is satisfied, then the compatibility conditions are 
satisfied for every smooth enough (C r ) desired output trajectory y<i{t). 

In addition, $*(t*+i,t*:) o Sk C Ck+i implies that 

Sk + 1 = $k{tk+i,tk) ° Sk 

is nonempty since 0 € Sq. Further dim(«S’„) = dim(£o,«)- 
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The additional condition for boundedness of solutions for t > Tj is met 
if and only if Sn. intersects £„,n- The linear variety Sn can be expressed as 

Sn = Itti(Sn) + v 

while the stable subspace of the internal dynamics can be expressed as 

£ a , n = /m(L S) Ar). 

It immediately follows from the previous considerations: 

Lemma 4 If assumptions 1-4 hold, then the exact output tracking problem 
with bounded solution is solvable if 

rank [ Sn Ln,» ] = n z (12) 

where n z is the dimension of the system internal dynamics at t = Tj. 

As a last remark to this section, note that if not only the internal dynamics 
subspace remains constant across the switchings as ensured by Assumption 
4, but also its stable and unstable subspaces do, then condition (12) is always 
satisfied and the problem has a solution for every admissible yd{t). Moreover, 
this solution is unique for every given yd(t). This implies that the exact- 
tracking problem for a system without switches is always solvable. 


3 . y d GIVEN THROUGH AN EXOSYSTEM 

In this section we consider how the previous results can be specified when 
the reference trajectory is not completely general but it is generated through 
a known linear exosystem, given by: 


X e — Ae(A:)x e 
Vd ~ Ce{k^jX e 


(13) 


First we solve the tracking problem when the initial state of the exosystem 
(at time Tf) is known in advance, hence the obtained result will be valid for 
the particular reference trajectory determined by the initial condition. Later, 
we will extend the approach to the case of unknown initial conditions. In 
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this case, we look for asymptotic tracking of the reference trajectories for 
arbitrary initial conditions of the exosystem. 

We begin by studying the case where the state of the exosystem z e (T t - ) is 
known. The system equation (4) becomes 

ii(t) = A y (k)x(t) + B,(k)C;(k)x.(t), (14) 

with A C*(k)x e (t) since all the time derivatives of the output can be 

written in terms of x e by using equation 13. 

The system evolution can then be rewritten as 

x{t) = T t )s{Ti) + H kt0 {t, Ti)x e (Ti ) (15) 

where 

k-i 

= h, k (t,tk) + k-l,j{tk—li 
i=« 

MM*) = /%^ )(t - T) B y (fc)C' e *(fc)4r fc ,o(r,7;)dr. 

Jt k 

and ^>k,o{x,Ti) is the evolution of the exosystem (analogous to equation 7). 
A solution to the exact tracking problem exists if lemma 2 is satisfied. The 
compatibility condition is satisfied for all initial conditions Xd(Ti) € Co if 
C* = C*. The stability condition becomes 

0 = siz ufi (Ti) + s 2 x e (Ti). (16) 

where 

[ 0 

■si = Z^N^f! NfliTfiT^Qa 1 0 

Iz Ut o(Ti) 

•S2 = Z U NHn,o(Tf,Ti) 

In what is to follow, we will assume that the above equation has a unique 
solution (iff sj is invertible). This yields a one to one relationship between 
the plant’s state and the exosystem as follows 

[ 0 

x d (Ti) = Qo 1 0 

. Zu,Ti(Ti ) _ 

0 

= Qo 1 0 AG(Ti)x e (Ti) 

. -Si 1 s 2 x e (T i ) _ 
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What is interesting is that we can also write the desired exact tracking 
state trajectory in terms of the exosystem state. Substituting the above 
expression in equation (15) we get 

x d (t) = G(t)x e (t) 


where 

am U)< 7 (B , + T,)c;] (17) 

It may be verified that G(t) satisfies the differential equation 

G k (t) = A,(k)G h (t ) - Gl(i)A,(k) + C" 

In the case of no switching a constant solution always exists for the above 
Lyapunov equation provided the eigenvalues of the exosystem A e are different 
from the zeros of the plant eigenvalues of A y . The above equation also 
provides a control strategy when the exosystem states are not known. We 
estimate the state of the exosystem as x e and regulate the trajectory x d = 
G(t)x e . The stability of such a controller is studied in the next section. 

4. Stabilization 

If the state of the exosystem x e is not known, then we could estimate it as x e 
with || x e (t) — x e (t ) llj < K e e aet || x e (0) — x e (0) || 2 . We use x d (t) = G(t)x e (t) 
as the estimated desired state trajectory, and stabilize this trajectory by 
using the control scheme shown in Figure 3. Note that the feedforward used 
(equation 3) is completely specified in terms of the exosystem’s state-estimate 
as follows 

u d = F k ,ix d (t) + F ka y d k \t) 

= F kA G(i)x e (t) + F k>2 C;x e (t) 

A F k x e (t). 

The state equations are of the form 

X = A k x + B k (FfcXe + K(x - G(t)x e )) . 

We require that the system in each interval is either stable or stabilizable. 
For simplicity, we assume that A(k) + B(k)K is Hurwitz for all k - the 
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Figure 3: Trajectory Tracking with Exosystem 
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arguments remain valid if the system is stabilized through any other feedback 
control scheme. . 

The desired trajectory satisfies 

Xd = AkXd + B k F k X e . 

Let e := x — Xd. Then the difference of the last two equations yields 

e = (A k + B k K)e + ( B k F k + B k KG(t)) (x e - x e ) 

The exponential stability of the error dynamics system follows from the fol- 
lowing lemma 

Lemma 5 Given the system, e = Ae + u(t) where K\e ait < || e~ At |( 2 < 
K 2 e a2t and || v(t) || 2 < K 3 e~ azt , .with Ki,K 2 ,K 3 positive , then e = 0 is an 
exponentially stable equilibrium point provided a 3 > a 2 > 0 

Proof: Using the variation of constants formula 

|| e~ At e(t) || 2 < || e(0) + / 0 * e~ AT v(r)dT || 2 

< II e (°) ll 2 + II Jo e~ AT v(r)dT || 2 
^ II e (°) ll 2 + / II e~ AT v{r) \\ 2 dr 

J 0 

< ||e(0)|| 2 + K 2 K 3 [‘e^--dT 

Jo 

< II e(0) || 2 + K 2 K 3 — 

«3 ~ 0:2 


provided a 3 > a 2 . Hence 

*MI e a "e(t) || 2 < || e~ At e(t) || 2 < || e(Q) || 2 + K 2 K 3 —^— 

«3 “ «2 

the rhs being a constant. Therefore, for all e > 0 there exists a positive 
constant K such that 


e(t) || 2 < Ke-^ 


16 




Figure 4: Example 

5. Example 

Consider the flexible structure, cantilevered at the base and free at the top, 
shown in figure 4. It is modeled (with Finite Element Method) with a single 
flexural element. The degrees of freedom are the translational motion at 
the base (xx), and the top (x 2 ), and the rotation at the top, x 3 . Input is a 
translational force at the base and the system output is x 2 . The structure is 
loaded with a mass m t , which is changed at several instances. The flexural 
element has the following properties: mass 420; length 1; elastic modulus 1; 
and cross sectional area moment of inertia 1. The objective is to maintain 
the top of the structure along a prescribed trajectory so as to facilitate the 
transfer of the load. The equations of motion can be described by 

Mi + Sx = Bu 

where _ 
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x = [xi x 2 xz\', and B = [1 0 0]'. 

In the standard form x = [x , x ]' (abuse of notation) we get the dynamics 
as 

x = AkX + Bkii 
Vk = C k x 

where 


A* = 


0 I ' 
-M~ l S 0 ’ 


;B k 


o 

M-'B 5 


Cfc = [0 1 0 0 0 0]; 


The desired trajectory is generated by an exosystem of the form x e = A e x e , 
where 


A e = 


0 10 0 

-10 0 0 
0 10 0 

0 0 10 


and the desired output is t/j = [0 0 0 l]x e & C e x e . We also switch the 
exosystem to A e = 0 4x4 at the initial and final times T,- = 0 and T) = 2 * pi. 
The first two states of the exosystem form an oscillator and the second state 
is then integrated twice to obtain the desired output. Note that (A e , C e ) is 
observable. Hence the exosystem states can be estimated. In our simulations 
we ensure that the output trajectories have a compact support [0, 27r] by 
choosing initial conditions of the exosystem of the form x e (Ti) = [0; *; 0; 0]. 

We also switch the mass m t on the structure (see figure 4) to take values 
m t = 0 Vf 6 [0,7r/2], m t = 100 Vt € (7r/2, and m t = 10 Vt 6 

(1.57r,27r], which denotes the jumps in the system. Note that C k remains 
constant through the switches and hence the compatibility conditions are 
always satisfied. As illustrated in Section 3, the map (j(0) : x e (Ti) — * Xd(Ti ) 
is given as 

-4.6715 -3.3239 0 0 ' 

0 0 0 1 

6.9406 2.6672 0 0 

-0.5935 -1.4710 0 0 

0 0 10 

2.6851 -0.8211 0 0 . 

As an example we simulate the forward dynamics with the initial condition 
for the exosystem as [0 1 0 0]'. The corresponding initial system state for 
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exact output tracking is z<f(T t ) = [—3.3239 0 2.6672 — 1.471 0 — .8211]'. 
The simulation results are shown in Figure 5, where the exact tracking state 
trajectory is shown - this desired state trajectory yields the desired output 
with an error of 10 -6 for a motion of 2 units. This error is believed to be 
due to the numerical integration schemes. Further the initial conditions are 
large and unrealistic. The initial conditions of the system are typically not 
the same as the initial conditions of the desired state trajectory - this re- 
sults in initial transient errors. If the system is stabilized then these errors 
decay exponentially, even if the system dynamics has switches. This ability 
to maintain tracking is a major advantage of our approach. Further, pre- 
actuation techniques to achieve these initial conditions (with output-error 
maintained at zero) has been developed in ( Devasia et al. (1996)) and we 
expect to integrate the two approaches in the future. 

Conclusions 

The problem of achieving exact output tracking for linear systems that 
present jumps in parameter’s values has been analyzed. We established nec- 
essary and sufficient conditions for the existence of exact output-tracking 
bounded state trajectories. When the reference trajectory is generated through 
an exosystem, the feedforward action needed to maintain exact tracking can 
be written as a time-varying feedback. Further, this time-varying feedback is 
related to a map from the state of the exosystem to the desired system state. 
The map is linear and is shown to satisfy an ordinary differential equation. 
For the case of systems without switches the theory presented reduces to 
the standard regulator theory. We also showed that the desired trajectory 
can be stabilized and presented the simulation results for an example flexible 
structure with switchingmass. 

Future work will attempt to remove the requirement of compact support 
for the output. There is also a need to address the tracking problem for 
systems whose internal-dynamics may not be hyperbolic. 
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Optimal Output Trajectory 
Redesign for Invertible Systems 


S. Devasia* 

University of Utah, Salt Lake City , Utah 84112 

I. Introduction 

G IVEN a desired output trajectory, inversion-based techniques 
find input-state trajectories required to exactly track the 
output. 1,2 These inversion-based techniques have been successfully 
applied to the endpoint tracking control of multijoint flexible ma- 
nipulators in Ref. 3 and to aircraft control in Ref. 4. The specified 
output trajectory uniquely determines the required input and state 
trajectories that are found through inversion. These input-state tra- 
jectories exactly track the desired output; however, they might not 
meet acceptable performance requirements. For example, during 
slewing maneuvers of flexible structures, the structural deforma- 
tions, which depend on the required state trajectories, may be unac- 
ceptably large. Further, the required inputs might cause actuator sat- 
uration during an exact tracking maneuver, for example, in the flight 
control of conventional takeoff and landing aircraft. 5 In such situ- 
ations, a compromise is desired between the tracking requirement 
and other goals such as reduction of internal vibrations and preven- 
tion of actuator saturation; the desired output trajectory needs to be 
redesigned. 

Here, we pose the trajectory redesign problem as an optimization 
of a general quadratic cost function and solve it in the context of 
linear systems. The solution is obtained as an off-line prefilter of the 
desired output trajectory. An advantage of our technique is that the 
prefilter is independent of the particular trajectory. The prefilter can 
therefore be precomputed, which is a major advantage over other 
optimization approaches (see Ref. 6 for further references). 
Previous works have addressed the issue of preshaping inputs 
to minimize residual and in-maneuver vibrations for flexible struc- 
tures; see, for example, Refs. 6 and 7. Since the command preshap- 
ring is computed offline, in Ref. 8, the use of noncausal prefilters has 
been suggested — such noncausality is allowable since the command 
preshaping is computed off-line. Further, minimization of optimal 
quadratic cost functions has also been previously used to preshape 
command inputs for disturbance rejection in Ref. 9. All of these ap- 
proaches are applicable when the inputs to the systems are known a 
priori. Typically, outputs (not inputs) are specified in tracking prob- 
lems, and hence the input trajectories have to be computed. The 
inputs to the system are, however, difficult to determine for non- 
minimum phase systems like flexible structures. One approach to 
solve this problem is to 1) choose a tracking controller (the de- 
sired output trajectory is now an input to the closed-loop system) 
and 2) redesign this input to the closed-loop system. Thus, we ef- 
fectively perform output redesign. 6 These redesigns are, however, 
dependent on the choice of the tracking controllers. 10 Thus, the 
controller optimization and trajectory redesign problems become 
coupled; this coupled optimization is still an open problem. In con- 
trast, we decouple the trajectory redesign problem from the choice 
of feedback-based tracking controller. It is noted that our approach 
remains valid when a particular tracking controller is chosen. In 
addition, the formulation of our problem not only allows for the 
minimization of residual vibrations as in available techniques 6 but 
also allows for the optimal reduction of vibrations during the ma- 
neuver, e.g., the altitude control of flexible spacecraft. 9 We begin by 
formulating the optimal output trajectory redesign problem and then 
solve it in the context of general linear systems. This theory is then 
applied to an example flexible structure, and simulation results are 
provided. 
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H. Problem Formulation 
System Inversion for Exact Tracking 
Consider a square system described by 

x = Ax + Bu\ y = Cx 

where x e 9ft", u e dt p , and y e SK 7 *. The inversion approach 2 finds 
a bounded input-state trajectory that satisfies the preceding system 
equations and yields the exact desired output, i.e., 

i ref = Ax rcf + Bu ff \ y d = Cx K f 

The inverse input-state trajectories can be described in terms of 
Fourier transforms as 1, 1 1 

Uff(J(0) = [COW - A)-' B]-'y d (jw) = G~ l (jco)y d {j co) 
XtMoS) = [ (J<oI - Ay l B]u ff (jw) = G x (ja>)u f f(Jt 0 ) 

This Fourier-based inversion approach has been extended to non- 
linear time- varying nonminimum phase systems in Ref. 2; however, 
we restrict our discussion to linear time-invariant systems. 

Remark We note two results. One, an inverse exists if the output 
and a certain number of its time derivatives are bounded. The number 
of time derivatives of the output that needs to be specified for an 
inverse to exist is well defined and depends on the relative degree of 
the system. 212 Second, for linear hyperbolic systems, if the inverse 
exists, then it is unique. 1,2 


Performance Criterion 

Trajectory redesign seeks a compromise between the goal of 
tracking the desired trajectory and other goals such as reducing 
the magnitude of input and vibrations. We formulate this redesign 
problem as the minimization of a quadratic cost function of the type 



u(t) T Ru(t) + x{t) T Q x x{t) 


+ [y(t) ~ y<im T Q,ly(.t) - *(i)]} d t 


where R, Q x , and Q y represent the weight on control input, states, 
and the error in output tracking, respectively. 

Using Parseval’s theorem, we rewrite our optimization problem 
in frequency domain as the minimization of the cost function 


/ oo 

{u{ja>y Ru{j(D) + x(ja>y Q x x(ja>) 

'OO 

+ ly(ju) - yd(J«>)YQyly(J<i>) - yaOV)]} d w (2) 

where the superscript * denotes complex conjugate transpose. 


Optimal Redesign of the Output 

Our main result is given by the following lemma, which shows that 
the optimal output trajectory redesign can be described as a prefilter, 
which maps desired output trajectories y d to the redesigned output 
trajectory y opl . This prefilter Gf does not depend on the particular 
choice of desired trajectory and hence can be precomputed. 

Lemma . The modified output trajectory y 0 pt is given by y 0 ptO^) = 
G f (jco)y d (jco), where 

Gf(jco) = 1 - G y [R + G X Q X G X + G;Q y G y ]~' 

x[r + g;q x g x ]g;' 

The modified input trajectory « opl is given by u m (Jw) = u ff(jco) + 
Vapiijco), where v„ pt (jw) = G v (jco)y d (ja>) and 

G v <jco) = -(r + g;q x g x + G;Q y c y y' 

x(r + g;q x g x )g (3) 

Note that the dependence on jco is not explicitly written for com- 
pactness. 

Proof. Without loss of generality, we rewrite the input u as the 
sum of the feedforward input, G~ [ y d , found from inversion of the 
desired trajectory, and an arbitrary v: 

u(jco) = Uff(jco) + v(jco) = G~\j(x))y d {j(o) + v(jco) (4) 

Substituting x(jco) = G x (jco)u(ja)) and y(jco) = G y (jco)u(jco) 
along with the preceding Eq. (4) for u into the cost function given 
by Eq. (2), we obtain 
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/ oo 

{[v + (R + G' x Q x G x +G* y Q y G y y X 

OO 

x (R + G' x Q x G x )G; l y d }* * (R + G;Q X G X + G\Q y G y ) 

x [d + (/? + G* Q X G X + GyQyGy'j 

x{R + G;Q x G x )G;'y d ] + (G; l y d y[(R + g;q x g x ) 

+ (r + c;Q x G x y * ( r + g;q x g x + c;QyG y y l 

x(R + G* x Q x G x )](G; l y d )}div 

Note that the cost function is quadratic in u. Therefore, the cost 
function is minimized by setting this quadratic term to zero, i.e., 
choosing v(jco) = v Q?t (ja>) = G v (jco)y d (ja >), where G v is defined 
by Eq. (3) in the lemma. The choice of v = v opt defines the optimal 
input M opt through Eq. (4) as 

WoptO'oO = [c;‘ (jco) + G v (jco)]y d (Ja>) (5) 

The result follows from y op i(j co) = G y (j co)u Qpl (j co) = [1 4- G y (j co) 
G v (ja))]y d {ja )). □ 


III. Example 

Consider a flexible structure consisting of two freely rotating 
disks connected by a thin shaft. A motor is attached between the 
connecting shaft and one of the disks. Input to the system is torque 
r provided by a dc motor, and the outputs are the angular rotations of 
the two disks 0i and 0 2 . These angular rotations are measured using 
potentiometers. The transfer function of an experimental system, 
which includes the rigid-body mode and one flexible mode, was 
obtained using a HP3562A Dynamic Signal Analyzer as 

0i 1.8139J 2 4- 0.3077^ + 6.1041 
t " j 4 + 0.2765.? 3 4- 6.1 04l5 2 
0 2 0.27s 2 - 0.1 187*4- 6.1041 

r ~ s 4 4- 0.2765s 3 4- 6.1041s 2 


with y s= 0 2 = [0 1 0 0]x. The control objective is to track the 
angular rotation 0 2 of the disk that is farthest away from the motor 
(see Fig. 1 for the desired output trajectory). 

The relative degree of a single-input/single-output linear system 
is the number of zeros at infinity. 12 For our system, with the torque 
as input and 0 2 as output, the transfer function has four poles and 
two finite zeros [see Eq. (6)]. Thus, the number of zeros at infinity 
are two, and hence the relative degree is two. This implies that the 
second derivative of the desired output, i.e., the desired angular 
acceleration profile of the output, uniquely determines the required 
input-state trajectoiy and the resulting structural vibration, 6\ — 02 . 2 
If the internal vibrations are to be reduced, then we have to relax the 
exact tracking requirement. Similarly, to reduce the required input 
amplitudes we have to compromise exact tracking. This tradeoff 
can be represented as the minimization of a general quadratic cost 
function (Sec. II) of the form 



u(t) T Ru(t) 4- x(t) T Q x x(t) 


+ [y(0 - yAt)] T Q y [y(t) - >,(/)]} dr 
where R = r, Q y = q y , and 

r i -i o o' 



0 0 0 0 


The scalars r, q x , and q y represent the relative weight on the re- 
duction of inputs, structural vibrations, and tracking errors, respec- 
tively. To reduce the vibrations and control inputs, we choose r = 1 , 
q x = 5000, and q y = 1 in our simulations. Figure 1 shows the modi- 
fication for a desired trajectory — about 10% of the final slew angle. 
The maximum magnitude of the required input, however, is reduced 
by 60%, and the corresponding structural vibration, 9\ — 02 , is re- 
duced by 20% (compared with results from exact tracking of the 
initial desired trajectory). 


With the state vector x chosen as x = [0! 0 2 9\ 0 2 ] r , the system 
equations can be represented in state-space form as x = Ax 4- Bu, 
i.e., 



ren 


‘ 0 

1 

0 

0 “ 

d 

@2 


-3.1555 

-0.1640 

3.1555 

0.3845 

dr 

0i 


0 

0 

0 

1 


-02- 


_ 2.8956 

-0.0899 

-2.8956 

-0.1124. 



r e, i 


" 0 “ 


02 

4- 

1.8139 

X 

0, 

0 


-02- 


0.27 


TV. Conclusion 

We formulated and solved the trajectory redesign problem in the 
context of linear invertible systems, including nonminimum phase 
systems. Thus, we provide a systematic approach to an optimal 
tradeoff between tracking desired trajectory and other goals such as 
vibration reduction and reduction of required inputs. The approach 
was applied to an example flexible structure, and simulation results 
were presented. Future work will address trajectory redesign for 
nonlinear systems. 
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Abstract 

In this paper we study the optimal redesign 
of output trajectory for linear invertible sys- 
tems. This is particularly important for track- 
ing control of flexible structures because the 
input-state trajectories that achieve the re- 
quired output may cause excessive vibrations 
in the structure. A trade-off is then required 
between tracking and vibrations reduction. We 
pose and solve this problem as the minimiza- 
tion of a quadratic cost function. The theory is 
developed and applied to the output tracking 
of a flexible structure and experimental results 
are presented. 

1 Introduction 

Large space structures, like manipulators used 
to assemble the space station, are lightweight, 
and hence flexible. The structural flexibility 
results in elastic vibrations. These vibrations 
are caused not only by exogenous perturba- 
tions but also arise during maneuvers like slew- 
ing. Hence the output tracking control for 
flexible structures is a formidable task. Re- 
cent works have solved the output tracking 
problem, for example, given a desired out- 
put trajectory, inversion-based techniques find 

*This work was supported by NASA-Ames Research 
Center Grant No. NAG 2-1042 


input-state trajectories that exactly track the 
output [1, 2, 3, 4]. These inversion-based 
techniques have been successfully applied to 
the control of multi-joint flexible manipulators 
in [5, 6], and to aircraft control in [7, 8, 9]. Note 
that the specified output trajectory uniquely 
determines the required input and state trajec- 
tories. Therefore, although these input-state 
trajectories exactly track the desired output, 
they might not meet performance requirements 
in flexible structures. For example, during 
slewing maneuvers of a flexible manipulator, 
the structural deformations - determined by 
the inverse state trajectories - may be unac- 
ceptably large. In such situations, a compro- 
mise is desired between the tracking require- 
ment and the other goals like the reduction of 
internal vibrations and prevention of actuator 
saturation - the output trajectory needs to be 
redesigned. 

Here, we pose the trajectory redesign problem 
as an optimization of a general quadratic cost 
function, and solve it in the context of lin- 
ear systems. The solution is obtained as an 
off-line prefilter of the desired output trajec- 
tory. Thus the redesigned output trajectory 
can be obtained through a simple convolution 
similar to the convolutions prevalent in the 
command shaping approaches, see for exam- 
ple [10, 11]. An advantage of our technique is 
that the pre-filter is independent of the partic- 



ular trajectory. The pre-filter can therefore be 
be pre-computed which is a major advantage 
over other optimization approaches (see [11] for 
references). 

Previous works have addressed the issue of pre- 
shaping inputs, rather than desired outputs, 
to minimize residual and in-maneuver vibra- 
tions for flexible structures - see for exam- 
ple [11, 12, 13, 14, 15, 16]. Note that in these 
problems, the command preshaping is com- 
puted off-line, and hence noncausal pre-filters 
can be used [10]. Further, minimization of op- 
timal quadratic cost functions have also been 
previously used to preshape command inputs 
for disturbance rejection in [17]. A common 
problem with all these approaches is that they 
are applicable when the inputs to the systems 
are kn own a priori. In tracking problems, how- 
ever, outputs are usually specified and not the 
inputs : input trajectories have to be computed 
from the desired outputs. For nonminimum 
phase systems (flexible structures with non- 
collocated inputs and outputs) the inputs to 
the system are difficult to determine. One ap- 
proach to solve this problem is to: (1) choose a 
tracking controller - the desired output trajec- 
tory is now an input to the closed loop system; 
and (2) redesign this input to the closed loop 
system. Thus, effectively, we perform output 
redesign [11]. These redesigns are, however, 
dependent on the choice of the tracking con- 
trollers [18]. Thus the controller optimization 
and trajectory redesign problems become cou- 
pled - this coupled optimization is still an open 
problem. In contrast, our optimal output tra- 
jectory redesign is independent of the particu- 
lar choice of tracking controller. We can, there- 
fore, decouple the trajectory redesign problem 
from the tracking-controller design. Note that 
if a particular closed loop controller is chosen, 
our approach is still valid. The formulation 
of our problem not only allows for the mini- 
mization of residual vibrations as in available 
techniques [11], but, additionally, allows for the 
optimal reduction of vibrations during the ma- 
neuver, which is required in maneuvers like the 
altitude control of flexible spacecraft [17]. 


We begin by formulating the optimal output 
trajectory redesign problem, and then solve it 
in the context of general linear systems. This 
theory is then applied to an example flexi- 
ble structure and experimental results are pro- 
vided. 


2 Problem Formulation 


System Inversion for Exact Tracking 

Consider a square system described by 
x = Ax + Bu 

y - Cx 

where x € 5R n , u £ $ft p and y € The inver- 
sion approach [4] finds a bounded input-state 
trajectory that satisfies the above system equa- 
tions, and yields the exact desired output, i.e. 


X re f — AXj.gf -]- Buff 
Ud ■— Cx re f . 


The inverse input-state trajectories can be de- 
scribed in terms of Fourier transforms as [1, 2] 


Uff(jw) = [C(jwl - A) 1 B] 1 yd{jw ) 
:= Gy 1 (jw)yd(jw) 


x ref(j w ) 


= [( jwl - A ) 1 B]uff(jw) 

■ = G x (jw)u f f(jw) . 


This Fourier based inversion approach has been 
extended to nonlinear time-varying nonmini- 
mum phase systems in [4], however, we restrict 
our discussion to linear time-invariant systems. 


Remark I We note two results. One, an in- 
verse exists if the output and a certain number 
of its time-derivatives are bounded. The num- 
ber of time derivatives of the output that needs 
to be specified for an inverse to exist is well- 
defined and depends on the relative degree of 
the system [19, 4]). Second, for linear systems, 
if the inverse exists then it is unique [1, 2, 3, 4]. 



The Performance Criterion 
Trajectory redesign seeks a compromise be- 
tween the goal of tracking the desired trajec- 
tory and other goals like reducing the magni- 
tude of input and vibrations. We formulate 
this redesign problem as the minimization of a 
quadratic cost function of the type 

1 - 0 o {u(t) T Ru{t)+x(t) T Q x x(t) + 

[y(t) - yd(t)] T Q y [y(t) - j/d(t)]} dt 

where R, Q x and Q y represents the weight on 
control input, states, and the error in output 
tracking respectively. 

Using Parseval’s theorem we rewrite our opti- 
mization problem in frequency domain as the 
minimization of the cost function 

/ OO 

{ u(jw)*Ru(jw ) + x(jw)*Q x x(jw)+ 

-OO 

[y{jw) - yd(jw)]*Q y [ y{jw ) - yd(jw)] }dw 

where the superscript * denotes complex con- 
jugate transpose. 

Optimal Redesign of the Output 

Our main result is given by the following 
lemma, which shows that the optimal output- 
trajectory redesign can be described as a pre- 
filter, which maps desired output trajectories, 
yd , to the redesigned output trajectory, y opt . 
This pre-filter, Gj, doesn’t depend on the par- 
ticular choice of desired trajectory and hence 
can be pre-computed. 

Lemma The modified output trajectory, 
y opt , is given by 

y op t(jw ) = G f (jw)y d (jw) 

where 

Gj(jw) = 1 — G y [R + Gx^QxGx^r 

G/QyGy ]- 1 [R + G X *Q X G X ] Gy - 1 . 

□ 


Remark II We point out two extreme cases. 
First, if the weight on the tracking error is zero, 
Q y = 0, but R is positive definite then we ob- 
tain v = — {Gy^ydj = —u/f. This implies 
that the input u opt = Uff + v = 0, i.e., the best 
strategy is not to track the desired trajectory 
at all. It is more optimal to remain at the equi- 
librium point where the cost is zero. Second, 
if the weight on the inputs and states are zero, 
i.e., R = 0 and Q x = 0 but with Q y positive 
definite then y op t = yd- This implies that exact 
tracking is optimal, and the cost is again zero. 

3 Example 

Consider an experimental flexible structure 
which consists of two discs connected by a thin 
shaft as shown in figure 1. These two discs can 
rotate freely. 



Figure 1 Schematic Experimental Setup 

The transfer function of the system, which 
includes the rigid body mode and one flexi- 
ble mode, was obtained experimentally using 
a HP3562A Dynamic Signal Analyzer. Input 
to the system is torque, r, provided by a DC 
motor and the outputs are the angular rota- 
tions of the two discs 61 , 82 - These angular 
rotations are measured using potentiometers, 
and the transfer functions are obtained as 

01 1.0684S' 2 + O.ISm + 3.5953 
7 “ S 4 + .2765 S 3 + 6.1041S' 2 

0 2 .lSOO-S 2 - 0.0699,5' + 3.5953 
7 ~ S 4 + .2765^ 3 + 6. 1041 S' 2 ( ’ 


Proof: See [20] 



With the state vector x chosen as 

' *1 ' 

$2 

x e x 
®2 

the system equations can be represented in 
state-space form as 

Ax + Bu 
O 2 — [0 10 0]x 


0 10 0 
- 3.2 - 0.16 3.2 0.38 

0 0 0 1 ’ 

2.9 - 0.09 - 2.9 - 0.11 _ 

0 ' 

1.81 
0 ' 

0.27 _ 

The control objective is to track the angular 
rotation, 62 of the disc, that is farthest away 
from the motor (see figure 1). The desired out- 
put trajectory and its time-derivatives as pre- 
scribed as shown in figure 2. 


x := 

V ■= 
where 

A := 
B := 


as input and 0 2 as output, the transfer func- 
tion has four poles and two finite zeros (see 
equation(3)). Thus the number of zeros at in- 
finity are two and hence the relative degree is 
two. This implies that the second derivative 
of the desired output, i.e. the desired angu- 
lar acceleration profile of the output, uniquely 
determines the required input-state trajectory, 
and the resulting structural vibration, 9\ — 0 2 
[4]. If the internal vibrations are to be reduced, 
then we have to relax the exact tracking re- 
quirement. Similarly to reduce the required 
input amplitudes we have to compromise ex- 
act tracking. This trade-off can be represented 
as the minimization of a a general quadratic 
cost function (section 2) of the form 


f-00 u(t) T Ru(t) + x(t) T Q x x(t)+ 

(y(t) - yd{t)) T Qy ( y(t ) - y d (t )) dt 


where R = r, Q y = q y , and 


Qx 


q x 


1-10 0 ' 

-1 10 0 

0 0 0 0 

0 0 0 0 



■ A \ 

1 1 — 

A 

■ 

■ ■ ■ desired velocity (degree/s) 


- 

—desired acceleration (degreed) 

1 1 « — 

, 


— upstieuapcaierayon tuey leery ^ , 1 , , j_ 
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time 


Figure 2 The Desired Output 

The relative degree of a single-input single- 
output linear system is the number of zeros at 
infinity [19]. For our system, with the torque 


The scalars r, q x , q y represent the relative 
weight on the reduction of inputs, structural 
vibrations and tracking errors respectively. In 
order to reduce the vibrations and control in- 
puts, we choose r = l,q x = 5000, and q y = 1 
in our simulations. Figure 3 shows the modifi- 
cation of the desired trajectory - about 10% of 
the final slew angle. 



Figure 3 Simulation: Re-Designed Output 







Figure 4 Optimal Input 



Figure 5 Simulation: Twisting of the Flexible 
Shaft 

The maximum magnitude of the required in- 
put, however, is reduced by 60% and the struc- 
tural vibrations, 0 \ — 62 , is reduced by 20% as 
shown in figures 4 and 5. The input-state tra- 
jectories found through inversion as described 
in equation 1 : two simulations are performed, 
first with the original desired trajectory, yd 
and second with the modified output trajec- 
tory, y 0 pt- These state trajectories are stabi- 
lized through feedback (see control scheme in 
figure 6), and preliminary experimental results 
are presented in figures 7 and 8. 



Figure 6 The Control Scheme 

With the reference feedforward inputs as the 
original inversion input Ufj and the modified 


(reduced) optimal input u opt (figure 4), the ex- 
perimental results verify that the small modi- 
fication of the output trajectory (figure 7) re- 
duces the maximum internal vibrations (50%, 
figure 8). Further reduction in the vibrations 
are possible if the tracking requirement ( Q y ) is 
relaxed. Also, it is possible to stabilize the sys- 
tem first through the feedback loop and then 
apply the trajectory redesign scheme to the 
closed loop system - but the redesign for our 
control scheme will not change since the input 
state trajectories always satisfy the state equa- 
tions and hence in the redesign the stabilizing 
feedback is zero. For example, («//, x r e /) satis- 
fies the state equations and thus the stabilizing 
feedback is zero if x = x re f and does not affect 
the output redesign. 




Figure 8 Experimental Results: The Internal 
Vibrations 

4 Conclusion 

We formulated and solved the trajectory re- 
design problem in the context of linear invert- 
ible systems - including nonminimum phase 








systems. Thus we provide a systematic ap- 
proach to an optimal tradeoff between track- 
ing desired trajectory and other goals like vi- 
bration reduction and reduction of required in- 
puts. The approach was applied to an example 
flexible structure and experimental results were 
presented. Future work will address trajectory 
redesign for nonlinear systems. 
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